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Abstract 

We have mapped the physics of polymer melts onto a time-dependent Landau- 
Ginzburg {tpl^ field theory using techniques of functional integration. Time 
in the theory is simply a label for the location of a given monomer along 
the extent of a flexible chain. With this model, one can show that the limit 
of infinitesimal concentration of a polymer melt corresponds to a dynamic 
critical phenomenon. The transition to the entangled state is also shown to 
be a critical point. For larger concentrations, when the role of fluctuations 
is reduced, a mean field approximation is justifiably employed to show the 
existence of tube-like structures reminiscent of Edwards' model. 
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Typeset using REVT^ 



Issues in polymer melt physics have continued to provide an enduring source of theo- 
retical investigations, owing to the complexity of the field. The basic physics of flexible 
polymer melts is embodied in Flory's theoremEl. This theorem states that for low densities, 
fluctuations in the melt are quite pronounced, while for large densities, fluctuations are so 
suppressed that the polymers behave as independent Gaussian chains again. One purpose 
of this paper is to modify the current perspective of polymer melt physics in the interesting 
low density regime. Rather than this state being close to a critical point in the static sensea, 
the connectivity of the chains in the system implies that the physics can be mapped onto 
a dynamic critical phenomenon. In the old picture, global properties such as the radius 
of gyration Rg (end-to-end distance) were shown to possess scaling properties, while local 
properties of chains could not be addressed. The new view-point allows a computation of 
the local structure of individual chains as well. Our estimate of the Flory exponent u is ap- 
proximately 0.631, in line with universality arguments, and differs from the usual estimate of 
0.588. More generally, the paper provides a powerful new method to study polymer physics. 
As an example we show that it can be used to describe the approach to entanglement as 
yet another critical point. This transition is beyond the scope of the standard tube model, 
a mean field approximation. 

The self-consistent field theory (SCF) discussed by de GennesS is a mean field approxi- 
mation to study melts of flexible polymers. Generalizing this model beyond its mean field 
roots has obvious advantages. Such a generalization has been attempted!. However, while 
this generalization was appropriate for the intended application, it ignored an important 
property of polymers. The theory involves a description of polymers in terms of a field '{/'(r), 
where r*is the location in physical space of any segment, such that |'?/'(f^)p is the probability 
of finding a segment at f. If the polymer is segments long, there is no representation in 
this theory of which segment (1 through N) this field refers to. In other words, reference to 
the connectivity of the chains is missing. This can be achieved using ideas from functional 
integration. 

The propagator for a single flexible chain may be represented byi: 
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where b is the bond length of the polymer, and where c}„ = This expression is 
obtained by considering only the entropy of a flexible chain. 

Alternatively, one knows from methods in functional integration thati: 



'Xn\ 




ij*{Ri, n)V'(i?2, 0) exp 



^^{x, n') 
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where 'D'^ip = 'D%l)*'Dil) , (3 = -j^, ks is Boltzmann's constant and T is the temperature. Thus 
we have another way of thinking about a system of flexible polymers, in terms of ip{x, n) and 
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an energy functional f3T which is isomorphic to one that describes diffusion. Here (x, n) 
labels the location x in physical space, of the n-th segment of a chain, and \%l){x,n)\^ is 
the probability of finding a polymer segment at a given location in space, as suggested by 
Eqn. ([I5|) below. In this sense we have a density functional theory in the style of Kohn and 
Sham. Kleinertji and Semenov et ali have proposed similar formalisms. However this paper 
uses the new formalism to probe the physics of polymer melts in a broader sense. 

It is possible to derive from the partition function Z = J'D'^ip exp— a 2p-point 
correlation function, which decouples at the non-interacting level into a product of Green's 
functions for p independent polymers. The main advantage of the functional path integral 
formalism is that one can now model more easily interactions in systems with large numbers 
of polymers. The following model, written to look like a Kohn-Sham type density functional 
theory describes excluded volume effects: 

Ho^ H = Ho + V - fi 




V=-mx,n)\' (3) 

where v is the usual excluded volume interaction parameteiH, and fi is the chemical 
potential invoked in the form of a Lagrange multiplier to conserve the number of polymer 
segments in the system. 

□ 

The self-avoiding walk of a solitary chain can be modeled by starting from Eqn. (0)0, by 
adding to the argument of the exponential on the right hand side, a series of terms which 
describe the excluded volume interaction between polymer segments. Caution must be used 
to apply this approach to a system of many polymers. 

To understand the physics in the \ip\^ model, let us extremize the functional density JF, 
while considering only solutions homogeneous in space and in n. This yields a maximum at 
ip = 0, and minima at: 

I^P = ^ (4) 

This of course leads to an infinity of solutions, equivalent within a phase difference. 
Following quantum field theory ideas, the phase of the field can be related to scattering, or 
interaction effects. The physical system can be thought of as localized regions where Eqn. 
is satisfied, separated by domain walls which permit the transition from one minimum to 
another. 

Note that Eqn.@ is equivalent to an estimate of the chemical potential (/io = cov) 
if the average number density cq is known. From this view-point, cq ^ represents a 
system of polymers which approaches a critical point from below. Here we have in mind an 
analogy with the usual Landau- Ginzburg 0^ model, where the vanishing of the coefficient 
of the quadratic term in the energy functional leads to a single well potential, signifying a 
critical point in the phase diagram. This issue was treated by de Gennes by mapping the 
polymer problem onto one in a zero-component (pj field theory (the self-avoiding walk of a 
solitary chain)i. The different perspective offered by our theory is that the physics of low 
concentration melts is really an issue in dynamic critical phenomena, given the degree of 
freedom represented by the variable n. 
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We can compute a dviiamic correction to the free particle Green's function, from the so- 
called Saturn diagramEniy. This diagram is the lowest order non-vanishing k, cu dependent 
contribution to the self-energy. To do this calculation, it is first convenient to write f3J^ in 
a dimensionless form, using ijj -i/^/i/cQ, and scaling all length scales by Cq^^^. This yields 
V —>■ a = cqv. In the limit of small fiQ, 

Goik,co)^ {G,\k,u)-^ik,uj)y' 

S fc,o;>«48a' / — ^ / jj. — 5 

where GQ{k,uj) is the Fourier transform of the Green's function introduced in Eqn.(0), 
and where we have performed the frequency integrations involved in the diagram using 
the method of residues, and used the mean field value for the chemical potential /xq = 
cqv 0. This integral can be found in Hohenberg and Halperini and they use Wilson's 
renormalization scheme to analyze the properties of this integral. We have been able to 
evaluate this multi-dimensional integral analytically. The integrals were performed using 
the identity 1/t = /q°° c/A exp(— At), and introducing a change of variables to the center- 
of-mass and relative co-ordinates of ki,k2. This allows a separation of variables to occur, 
permitting an integration over the momentum variables. The final integration is performed 
using the identity exp{—Ax)x''dx = A^~^''r{l + k). In using this identity, we have to 
perform the integration first in arbitrary dimension d, where the identity is valid, and then 
continue it formally to d = 3 — e, e— >0. The infinities in the integral are then isolated into 
the Gamma function. As usual, the infinite part is removed by introducing an appropriate 
counter term in /?JF. This amounts to a renormalization of b. The finite result is: 

S(fc,cj)= a'^A{2mf~^^^ (-icu + k'^/{6m) - 3/io)^ (C - \n[2m{-iu; + k'^/{6m) - 3/io)]) 
C= ln[4V37r] + diGamma(3) ^ 4 

u)^ a^AC{2mf~^l^ (-zcj + A;V(M - 3/io)^'^^^ (6) 
where A = 2^/3^l~^^'^, and m = ^ 'fs 2 • '^^^ ^^^^ approximation in Eqn.(^) applies when the 

C(j b 

argument of the logarithm has a magnitude less than one. Since /xo ^ 1? the approximate 
scaling form thus holds for \uj\ < 1 and k < \^18/b, with an error of a few percent or 
less. In writing down these equations, we have implicitly performed mass renormalization!. 
The effect of the self-energy S is more pronounced at small length scales than it is at long 
wavelengths. The renormalized Green's function can also be written for — > as: 

-tuj + [-^Q^) k^ - /io) 



where 



jl{uj,l2o)= /io - a^A{2mf-^/^{-iiu - 3/io)^"^^' 
52(cu,/io)=6'co/' + A62(cu,/io) 

Ab^{uj,i2o)= -a\6A){2mf'^/^ (2 - 1/C) {-luj - 3/io)'"'^^ (8) 



4 



Note that this approximation holds in the regime discussed below Eqn.(^). Now —iuj = 
d/dn ~ — q\, with p, q referring to the p*^ and g*'* segments respectively, while k is the 
inverse separation in physical space of these segments. In the vicinity of the critical point 
/io = 0, displays a scaling property viz., ~ —\p — q\~'^ in the appropriate regime, where 
we have defined a scaling exponent a 

a = (1 - 1/C) ~ 0.75 (9) 

In this regime, the coefficient of the scaling term behaves as ~ a'''^, a weak dependence, if 
we estimate f ~ 6^. 

It follows from Eqn.(^ that the effect of the gradient-smoothing term in the energy 
functional is reduced by an amount — |A6^|, and this effect is dominant for segments sepa- 
rated by a relatively short distance, as it vanishes in the limit of infinite separation. It means 
that if segments on a chain happen to be in close proximity in the melt, they will tend to 
stay together due to the reduction of the smoothing term. Neutron scattering experiments 
or numerical simulations may provide verification of this notion. Neutron experiments are 
customarily employed in the long wavelength regime to investigate quantities like the radius 
of gyration of polymeric systems. Other techniques, such as those used by Smith et al0 may 
be more useful. 

The chain takes on the appearance of pearls on a string which push the ends of the 
string further away from each other. One way to calculate the scaling properties of the 
radius of gyration within the current model is to compute vertex corrections within the 
current model, and then implement Wilson's scaling arguments. However, it is easier to go 
back to the original theory defined by the energy density ip*{dn + H)iIj, and truncate it with 
dn —>■ 1/N where N is the average chain length of the polymeric system. This is justified on 
the grounds that we are interested only in long range fiuctuations. Accordingly, we also have 
to restrict the n' integration in /3JF to a small neighborhood around A^. After performing 
this renormalization, the entire machinery of the static theory of critical phenomena applies. 
And it follows that very near the critical point /iq = O"*", the correlation length ^ = i?^ ~ A^'', 
where the lowest order field field theoretic techniques yield u ~ 0.6, the value obtained by 
Flory. 

More accurate calculations yield! a value for z/ closer to 0.63lii. The standard model 
for examining the effects of excluded volume interactions yieldsB u ~ 0.588. The difference 
between this model and ours was discussed below Eqn.(^. 

There is a transition to an entangled state at a value of the concentration cq < 1/v. To 
see this let us begin by computing vertex corrections, which can be thought of as giving 
rise to an effective coupling constant. The lowest order correction comes from the so-called 
fish diagram with two internal lines. The frequency integral involved in the calculation is 
performed by closing the contour in the lower half-plane. The remaining momentum integral 
contains a divergent part which involves an integral from Y^2m/io to oo. This is written as 
an integral from to oo minus a finite part from to v^2m/io. The infinite part is removed 
in the usual manner with a counter term, and amounts to a renormalization of the coupling 
constant. This leads to a new coupling constant a^. 

aR{q)= a - a^f{q) 
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where q = (g, cUg), etc. and /iq << 1. Notice that to this order in perturbation theory, 
an increasing concentration signified by an increasing a leads to a lower effective coupling 
constant, consistent with the latter half of Flory's theorem. We can define a beta function 
as is done conventionally in Renormalization Group theory, viz.. 



d&R 



-Sir + &^ 



dhiy 

y= -iuq - 2/io = l/N - 2fio 



(6/ic 



N3/2 



aR{q = 0) 



(11) 



where y plays the role of a scaling parameter, and as y ^ Q {N ^ (2/Uo) we see that 
dtR fiows towards a nontrivial fixed point viz., otR =1. As an example, it can be checked 



in the case of polydimethylsiloxane (PDMS), that Cq ~ 3 x 10~^ cm~^, h ~ l.SA, yielding 
Nc ~ 167, which is close to the experimentally determinedEl critical entanglement chain 
length of ~ 200. Similarly, for polystyrene, using h ~ lA, we get Nc ~ 320, reasonably 
close to the critical chain length estimated by the viscosity measurements of Onogi et alcj. 
To the lowest order, the structure factor S{q) = (2/3)r(g). It diverges as at y = 0, 
analogous to the behavior at a critical point. The Young's modulus Y of the system is given 
byii Y oc S{q = 0,uj), and it will show a dramatic rise as the chain length N is increased 
towards for /iq 7^ 0. For these reasons we can identify this critical point with the onset 
of entanglement. Note that the approximation to S{q) we have used here will not work 
past the transition to entanglement. There is some indication of fiuctuating behavior near 
the entanglement transition in the classic viscoelasticity measurements of Onogi et a 10, as 
indicated by earlier phenomenological analyses^'S. 

Fluctuations in the system tend to decrease as Cq increases away from the critical point 
just discussed. So for finite concentrations we expect mean field theory to hold. In this 
limit, we get a time-dependent Landau- Ginzburg equation (time = -n.), by extremizing f3J^ 
with respect to ip*: 




ip{x, n) = 



(12) 



For the moment, let us suppose that there is no n dependence in ip. Then the model 
reduces to the SCF eouation. Noting the similarity of this equation to the one used to 
describe bosonic fiuidalZl, we seek tube-like solutions in cylindrical geometry, viz., ^/^(p, 0) = 
exp(is0)/(/9), s is an integer: 



pdp \ dp I \ 



/(p)-f(p) = o 



where p has been made dimensionless by the length scale 



a = b/ i/Scof 



(13) 



(14) 



We see that for p — > 0, assuming |/(p)| — ^ 0, we get a solution proportional to Js(p) ~ 
By choosing the amplitude of this solution correctly, one can match it for large p to 
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a constant solution /(p) = ±y 1 — s^/p^ (see Fig.(|I])). This solution may also be continued 
mathematically to s — 0. But then the s = solution has an infinite slope at the origin. The 
radial extent of these solutions is ~ s a, for s 7^ 0. From excluded volume considerations, the 
maximum density possible in the system is 1/f , and for densities approaching 1/f , a ~ 6/ a/6. 
Physically, these solutions are reminiscent of the tube model of Edwards^. They indicate 
that the tubes are not empty, but have a concentration of polymers given by [/(p)]^- It may 
be possible to use numerical simulations to check if the qualitative features of the occupation 
profile [/(p)]^ can be reproduced^. 

As discussed by Fetter, the higher the value of s, the higher the energy of the tube- 
like configuration, and the system will prefer to have many tubes with a low s than a few 
with a large s. We have thus derived from our theory a tube model with a tube radius a, 
which is exceedingly large for small densities, and asymptotes smoothly in an inverse square 
root fashion to ~ as the density approaches 1/v. Since the monomer length b is of 

the order of a few Angstroms, the magnitude of the tube radius for cqV ~ 10~^ is a few 
monomer lengths. This is in qualitative agreement with estimates which can be found in 
the literature^. 

In the theory of superfiuidslll, these cylindrical structures are viewed as idealizations 
of vortices, caused by the rotation of the fiuid. In our case we can derive an equation of 
continuity from Eqn.([T^: 



dm 



2 



dn 

r- 



-V-3 + S 



(^*VV' + c.c. 



S= (^^)|V^I' + 2po|^P-^|^r (15) 

It follows by examining the current j, that the tube-like solutions are not caused by 
rotation, but rather, there is a radial velocity due to the p dependence of the solution. There 
is yet another solution to Eqn. ([T^) , and is obtained by assuming that it is dependent solely 
on n, which yields ip{n) = [1 + exp(— 2na)]~^/^. More general solutions can undoubtedly be 
found numerically by studying Eqn.([T^) with various boundary conditions. 

We speculate that the phase of the field ip may be useful in quantifying the notion of 
entanglement in polymers. Ultimately the theory needs to be generalized to handle polymer 
dynamics. 

This research is supported by the Department of Energy contract W-7405-ENG-36, under 
the aegis of the Los Alamos National Laboratory LDRD polymer aging DR program. 
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FIGURES 

FIG. 1. Numerical solutions of Eqn.(|l3[) for s=l and s=2. Both curves asymptote to unity, but 
on different scales. 
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Chitanvis, Fig. |l] 
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